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Abstract: 

We calculate radiative corrections to the light-cone sum rule for the semileptonic form 
factor in B —>■ -nev decays and thus remove the major uncertainty in determining the CKM 
mixing angle \ V u b\ by this method. We discuss the remaining uncertainties and perspectives 
for further studies. The structure of the radiative corrections suggests factorization of soft 
(end-point) and hard rescattering contributions in heavy-to-light decays in the heavy quark 
limit. 
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1. Precise determination of the CKM mixing angle \ V u b\ will most likely come from semilep- 
tonic B decays. Two competing strategies exist, the study of the end-point region in in- 
clusive decays and the analysis of various exclusive decay channels, most notably B — > neu 
which is the easiest one to access experimentally. The theoretical challenge in the latter 
case is to calculate the B — ► 7r transition form factor f+{q 2 ) induced by the weak vector 
current: 

(7r(p w )|«7 M 6|S(p B )> = (PB+Pn)»f+((1 2 ) + (PB ~Pn)^f~(q 2 ) • (1) 

Here q 2 = (j?g — p n ) 2 is the square of the momentum transferred to the lepton pair; the 
second form factor /_ does not contribute to the decay rate for zero mass leptons and 
will not be considered in this letter. The decay rate is dominated by contributions from 
small values of q 2 , for which the u quark produced in the weak decay has large energy 
of order m b j2 in the B meson rest frame (m& is the b quark pole mass). A consistent 
QCD description of such processes exists up to now only in the Sudakov limit in 
which case contributions from large transverse quark-antiquark separations are suppressed 
and the form factor is dominated by hard gluon exchange ||. This limit is theoretically 
interesting, but not relevant for realistic b quark masses. In fact, we will argue that "soft" 
contributions related to large transverse distances exceed the "hard" contributions by an 
order of magnitude (and have opposite sign). Any theoretical approach to heavy-to-light 
decay form factors aiming at quantitative predictions has to deal with these "soft" terms 
explicitly, as do e.g. light-cone sum rules and various models ||. 

The consistent separation of "hard" and "soft" contributions in heavy-to-light decays 
presents an unsolved problem which, to our knowledge, has never been tackled in a sys- 
tematic wayQ. In this letter we examine a possibility for such a separation, based on our 
calculation of the radiative corrections to the corresponding light-cone sum rules. 

The basic idea of the light-cone sum rule approach is to consider a two-point correlation 
function replacing the B meson state by a suitable interpolating current: 

n>|,g 2 ) =i I rf 4 xe-^(vr(^)|T{M(0) 7 ^(0)&(x)^5^)}|0) 



(p fl + p.) M 11+ (p|, q 2 ) + (p B - U_{p 2 B , q 2 ). (2) 



Tl + (PW) = fBml J 2 +iq2 \ +-, (3) 



The Lorentz- invariant function 11+ has a pole at p 2 B = m 2 B corresponding to the contribution 
of the B meson: 

m b m B - p B 

where me is the B meson (pole) mass and the dots stand for contributions from higher- 
mass resonances and the continuum. The B meson decay constant fg is usually defined 
as 

m b (0\di l5 b\B) = m 2 B f B . (4) 

On the other hand, the correlation functions in (0) can be calculated in the Euclidean 
region, p\ — m 2 B large and negative, using the light-cone expansion. Up to higher twist 



1 One reason being that in the classical application, the pion form factor, the "soft" terms are suppressed 
by one power of the momentum transfer. 
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corrections the product of b quark fields can be substituted by a perturbative propagator 
and the tree-level result is 



n+(p|, q 2 ) = \Urn h ^ du 2 M Y ] _ 2 . (5) 
2 Jo mi — up B — uq 2 



Here and below u = 1 — u, (j) n (u, p) is the leading twist pion distribution amplitude defined 
by 

(7r(p 7r )|w(0) 7i , 75 d(x)|0) x =° -if vPm> f'due^Mu^) (6) 



o 

in Fock-Schwinger gauge and p is the factorization (renormalization) scale. The variable u 
has the physical meaning of the momentum fraction carried by the u quark in the infinite 
momentum frame. 

Next, we use the concept of duality, assuming that the contribution of the B meson 
to (|3|) corresponds to an integral over the spectral density calculated within the light-cone 
expansion, Eq. (§), in a certain duality interval: 

n?»-» (pW)= r^ P (s, q 2 ). (7) 

Jm z b S — p B 

The width of the duality interval is characterized by the parameter sq (continuum thresh- 
old) and in general is of order s ~ m \ ~ 0(m^). The spectral density p(s, q 2 ) immediately 
follows from Eq. (|5]) after a simple change of variables s — > (m 2 — q 2 )/u + q 2 . Thus, in this 
particular case the duality restriction on the maximum invariant energy so translates into 
a restriction on the minimum momentum fraction Uq carried by the u quark: 

m 2 b -q 2 

u = 5"- (8) 

so - q 2 

Equating the representations @ and (|5|) and isolating the contribution of the B meson 
according to (|7|), we obtain the simplest version of the light-cone sum rule for f + (q 2 ), which 
neglects higher twist and radiative corrections. In order to suppress the contributions of 
higher order states, it is written in a Borel transformed form, which amounts to replacing 
the factor l/(s — p 2 B ) by exp(— s/M 2 ) and l/(m 2 B — p 2 B ) by exp(— m 2 B /M 2 ): 

fsm 2 B Jrn 2_ m 2 n)/M 2 , 1 r 1 6 v {u,y) 



(m 2 b -m%)/M s = t / du I^irl e -«K-g 2 )/( U iVf 2 )_ ^ 



f n ml ' ■>-<-■*< 2 ./, (: u 

Here M 2 is the Borel parameter. 

2. The accuracy of the sum rule (|Sp can be improved by including higher twist and 
radiative corrections. The former ones have been calculated earlier 0, [J; Eq. (79) in 
gives the complete result to twist 4 accuracy, which we will use in the numerical analysis. 
We have calculated first order radiative corrections to the leading twist term ([|) The 
calculation is straightforward and similar to earlier calculations of the radiative correction 
to the 7* — ► 7T7 form factor J7|, but more cumbersome because of the nonzero b quark 
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mass. In the MS-scheme, the result reads: 



fB m B c (m 2 b -m%)/M 2 



/+ (g 2 ) = - td J^^K -^i- 

2 JUQ u 



-<j 2 )/(uM 2 ) 



C F ^ ^ { I" du MU ^ ] r dt ^4 (- + Li) a* ' 

4-7T I Jo U Jm\ t - q 2 \t 



ut-\-uq -m. 



du 



u 



«0 



— / dt(mi-t) 



mf — ut — uq 1 



r 1 fafan) _ziz^! 

+ du - e «Af 2 

Juq U 



ut-\-uq —Tn, 



SO 



dt {m 2 b - t) 



m{ — ut — uq 2 



u 



so 



dt U 



*° , .ml — q 2 



dt- 



t-q 2 



m 2 h — t 
t(ml — q 1 



r 1 Mu,fi) _z^l 

+ du e u*f a 



5 j E uM 2 3 wM 2 I fm 2 -i 
— + 2 In — ^ -in — — + -Ei 



mi 



uM 2 



r s roc f m 2 _ „2\ T 

+ dtL 2 + dtX^ — + 



/•to 

+ / dt 



ml — ut — uq 2 
7 > + (ml - q 2 ) [L x + L 2 ) 



*0 /l \ 

(1/ I - — L 2 I o 



1 (t - m 2y2 

2 2^2 



e um 2 — 1 
m 2 , — t 



00 dt 



/„ mf-t 



(t ~ m 2 ? 
2t 2 



(ml-q 2 ) (L 1 + L 2 ) 



(10) 



where mj is the one-loop pole mass, Ei(x) the exponential-integral function, defined as 
Ei(a;) = — J^° x dy e~ y /y, t = usq + uq 2 and 



t-q 2 



tfi 2 



-1+ln 



t-q 2 



mi 



2\2 



+ ln 



(t-m 2 ) 
t/i 2 



(11) 



Details of the calculation will be published elsewhere. We have checked that the /i de- 
pendent terms cancel the \i dependence of the distribution amplitude 4> n (u, /i) to leading 
logarithmic accuracy ||. We have also checked that our results agree with the recent re- 
sults of Ref . |J , where a similar calculation is reported, both analytically for the amplitude 
in momentum space and numerically for the final sum rule after Borel transformation and 
continuum subtraction. Our representation for the final answer is simpler. 
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3. In order to gain a better understanding of the structure of radiative corrections, it 
is convenient to go to the heavy quark limit m& 
parameters fs, sq and M 2 have to be rescaled as 



oo. In this limit the dimensionful 



/ B = l/ v ^r at (/i = m b ) 



so 



ml + 2m b uJo, 



rris — rrib = A, 
M 2 = 2m b r, 



(12) 



where uq and r are the nonrelativistic continuum threshold and the Borel parameter, 
respectively. The heavy quark limit of the form factor depends crucially on the value of 
q 2 . As we are only interested in small values q 2 <C m 2 , we will set q 2 = in this discussion 
for simplicity. 

In this limit the momentum fraction cut-off @ becomes Uq = 1 — 2uJo/m b , so that the 
integration region in (|3|) shrinks to a narrow interval near the end-point, corresponding 
to all the pion momentum being carried by the u quark. Since close to the end-point 
</> n (u, fi) ~ (1 — u), we obtain 



-A/r 



mbfTT 



/+(0) 



nn, 



(1) / duJUJQ-^^ 



-3/2 



(13) 



at least 



where </4(u) = (d/du)^> n (u). Note that </4(l) < 0. One thus finds /+(0) ~ m b 
at tree level [g 0]. 

The full expression for the radiative corrections looks rather complicated; to make its 



structure more transparent we take the limit r 
duality approximation: 

/ stat (m fe ) r _ 3/2 



oo, corresponding to the so-called local 



\m 



/+(0)] 



-^o 2 C(i; 



71 



h In 

t 2Wn 



l ln2 m^ l ln 2^ 
2 2u 2 y, 



7T 



1 - In 




2u; r 



-ln^ I'du^- . (14) 

fl J JO \ U ) [I Jo u 

This expression deserves to be studied in some detail. Let us interpret the two pieces: the 
first term on the right-hand side must be identified with the soft (end-point) contribution 
including the Born-term and its radiative correction, while the second term corresponds to 
the usual mechanism of hard gluon exchange. 

The dependence on the collinear factorization scale \x must cancel the scale dependence 
of the pion distribution amplitude. This implies that the structure of terms in ln/i in 
the hard contribution is fixed by the structure of the leading order soft term which is 
proportional to (p' n (l,fi). Indeed, we find 

d a d r /"i 

7i — Ki 1 ^) = — c f~r / d v v °( x > v) &r(f> v) 

dm /j, n ax Uo 




V(n) + W0^(1) (j) n {u) 



v? 



u 



(15) 
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where V (x, y) is the usual Brodsky-Lepage kernel, so that the structure of In \i terms in ( |T4"D 
is reproduced. Note the subtraction term accompanying the naively divergent expression 
/ du (j) n (u, fi)/u 2 , which is similar to the usual "plus" prescription in the evolution kernel. 

Local duality means that we identify the B meson with a h quark accompanied by 
an arbitrary number of light quarks and gluons with total energy less than ujq (in the b 
quark rest frame). Consider the "deep inelastic" cross section of neutrino scattering off 
a pion, in which one selects the contribution of the charged weak current with a b quark 
in the final state: da j dM\{y e + 71 — > X b + e), where M\ is the invariant mass of the 
hadronic final state. In the approximation adopted in this letter we identify the integral 
of this cross section over the region of small invariant masses M\ < m 2 + Im^io® with 
the square of the form factor for the inverse process Fivix — > eB) (up to kinematical 
factors). This interpretation is useful in several aspects. It is easy to check that the "soft" 
contribution in ([L4]) corresponds to the would-be leading twist contribution to the deep 
inelastic cross section, while the "hard" contribution (involving interaction with the quark- 
spectator) is of higher twist. These two terms are of the same order (in the b quark mass) 
since the leading twist contribution is additionally suppressed by a factor 1/m 2 for small 
values of the invariant mass of the hadronic system, because the corresponding parton 
distribution (in the pion, for the case at hand) vanishes at the end-point x — ► 1. With a 
low collinear factorization scale \i ~ u>o all quark mass dependence in ( |I4"D is due to the 
"soft" contribution and thus to the leading twist part of the cross section. This suggests 
that resummation of heavy quark mass logarithms can be done using the same techniques 
as for the end-point spectrum in inclusive b — > u decays [[II]. A detailed discussion goes 
beyond the scope of this letter. 

A final remark concerns the size of the radiative correction. With the natural factoriza- 
tion scale \x = 2u ~ 2 GeV and with mj ~ 5 GeV the quark mass logarithms are of order 
unity and the large constant term dominates. It has to be compared, however, to the large 
radiative correction which was found in QCD sum rules for the decay constant J stat . In 
the same (local duality) approximation and neglecting contributions of condensates, one 
finds HIT 



pstat ' UJ < 



n\m b ) 



' 2 \ . ®s n /15 1 2 3 m b 

1H Cfl — 71 H — In 

7T * V 8 6 4 2u 



71 



(16) 



We see that the large radiative corrections almost cancel each other between / stat and the 



right-hand side of (14). This shows that the form factor itself is free from large radiative 
corrections in the light-cone sum rule approach. The same cancellation takes place in the 
complete expressions with finite b quark mass. 

4. We have carried out a detailed numerical analysis of the complete sum rules with 
finite b quark mass and including radiative and higher twist corrections. To this end we 
substitute the value of Jb in fllPf ) by the corresponding sum rule [|ll| including radiative 



corrections with the same value of the continuum threshold and the same Borel parameter. 
In particular, we use the range 5 GeV 2 < M 2 < 8 GeV 2 and m;, = {4.6, 4.7, 4.8} GeV with 
the continuum thresholds Sq = {35, 34, 33.5} GeV 2 , respectively. The resulting value for 



fs is (175 ± 25) MeV, which is in agreement with the current lattice average fl4| . We do 



not attempt a renormalization group improvement (resummation of lnm&) of the sum rules 



and set the factorization scale to ji — Jm B — m b ~ 2.4 GeV. The scale dependence of our 
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results is in fact negligible. A remark is in order about the actual choice of M 2 in Eq. (|TO|) . 
The expansion turns out to be essentially in inverse powers of uM 2 rather than M 2 itself. 
In order to avoid small values of this "effective" parameter uM 2 for large q 2 (where u can 
become small), we use the rescaled value M 2 = M$ pt /(u) in the Borelized expression for 
the correlation function Il + and M| L in the above mentioned range in the sum rule for 
(u) is the average value of u in the integral in with (u) ~ 0.87 for q 2 = and (u) « 0.72 
for q 2 = 15GeV 2 . Variations of the Borel parameter in the chosen window result in an 
minimal uncertainty in the prediction for the form factor of about 10%, independent of q 2 . 
We are going to discuss the other remaining uncertainties in detail. 

Radiative corrections. Radiative corrections to the light-cone sum rule have been pre- 
viously suspected to be significant. It was expected, however, that large corrections to the 
correlation function n + partially cancel with large corrections to the coupling fg in the 
ratio /+ ~ II + // B . To quantify this effect, we write, schematically, 

n° ( 1 + ni 

f + „ _tA tW (17) 

where and /g are the tree-level contributions and and / B specify the corrections. 
For central values of the input parameters, we find fg ~ 3.0 and = {2.4,2.3,2.2,2.1} 
for momentum transfers q 2 = {0, 5, 10, 15} GeV 2 , respectively. Thus the corrections indeed 
cancel each other to a large extent. The "hard" contribution to 11+, defined by the terms 
in ( |T0D involving an integral over u in the interval < u < Uq, is —0.5 (at q 2 = 0), whereas 
the radiative correction to the "soft" contribution u < u < 1 is +2.9. The resulting net 
effect of radiative corrections is shown in Fig. [I], where we plot /+ as a function of q 2 
for two different choices of the pion distribution amplitude (DA), see below. The curves 
marked LO (NLO) are obtained by neglecting (including) radiative corrections in both the 
numerator and the denominator in (|i~?P . The size of the correction is at most -7%, thus 
providing an a posteriori justification of the procedure of Refs. f|, 0, || to use a low value 
of fs in leading order light-cone sum rules. 

Pion distribution amplitude. The main input in the sum rules is the pion distribution 
amplitude (p n defined in (||), which can conveniently be expanded in a series of Gegenbauer 
polynomials with multiplicatively renormalizable coefficients a n (to leading logarithmic 
accuracy) :Q 

fa(u, y) = Qu(l -u)[l + a 2 ( y u)C 2 3/2 (2u - 1) + a 4 (/i)C 4 3/2 (2u - 1) + . . .] . (18) 

The values of the nonperturbative coefficients a n can be restricted from experimental data, 
most notably from the 77*7r form factor. As argued in jT5], the recent measurements by 



CLEO [[16] are consistent with the distribution being close to its asymptotic form. On the 



other hand, the old QCD sum rule results |L7| indicate sizable corrections which, in view 



of the criticism raised in [fpjj], are probably overestimated. In our analysis we use the QCD 



2 We neglect two- loop anomalous dimensions and also the mixing of the a n , which occurs at two-loop 
accuracy. The corresponding expressions are available in the literature and can easily be incorporated in 
the analysis. Their effect is, however, negligible compared to uncertainties in the numerical values of the 
a n . 
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Figure 1: f+(q 2 ) as a function of q 2 for two different sets of the leading twist pion distri- 
bution amplitude. The effect of including radiative corrections is a reduction of the form 
factor by about (4-7)%. 
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Figure 2: The several contributions to the light-cone sum rule for the form factor as a 
functions of q 2 , using the leading twist distribution amplitude of |19 . 
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Figure 3: The spectrum dB(B — > 7ieu)/dq 2 as a function of q 2 . Solid line: central values of 
the input parameters, using the distribution amplitude (DA) |0|; dashed-dotted line: the 
same using the asymptotic DA; dotted line: the vector dominance approximation (VDM). 
The dashed lines indicate the range of the theoretical uncertainties. See also text. 



sum rule estimates a 2 = 0.35 and 04 = 0.18 at the scale fi = 2.4 GeV |]T9, g as upper bounds 
for possible corrections to the asymptotic distributions. The resulting q 2 dependence of 
the form factor is displayed in Fig. 1 for the two choices — asymptotic and QCD sum rule 
motivated — of the pion distribution amplitude. In Fig. 2 we show contributions of the 
asymptotic distribution and the corrections separately. It is seen that (a) the contribution 
of the second Gegenbauer polynomial ~ a 2 can reach 20% and significantly affects the 
shape of the form factor as a function of q 2 and that (b) with current estimates of their 
magnitude, the contributions of higher polynomials (~ a 4 ) are negligible. We stress that 
the size of these corrections will eventually be determined from experimental data. 

Higher twist corrections. From Fig. 2 it is apparent that the twist three contribution 
to the sum rule is large and of the same order as the leading twist contribution. An 
inspection shows that this large contribution is entirely due to the asymptotic two-particle 
distribution amplitudes of twist three, whose normalization is fixed by the quark condensate 
p0fl . A variation of (qq) within the conservative limits —(230 — 250 MeV) 3 yields an 
uncertainty in /+(0) of at most ±0.012, i.e. 4%. Corrections to the asymptotic two-particle 
twist 3 distribution amplitudes are related to the three-particle distribution of twist 3 and 
proportional to the coupling defined by the matrix element 

(0\da^ 5 G a u u\7r + (p)) = 2if 3wP ^ Pu + ... (19) 

The constant was estimated from QCD sum rules to be = 0.0026 GeV 2 ±30% (at 
the scale 2.4 GeV) JT7J. The corresponding contribution to the form factor is within 2%. 
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Finally, the twist 4 contributions also turn out to be unimportant numerically, see Fig. 2, so 
that we use the full set as specified in || without detailed error analysis. The only essential 
uncertainty in the higher twist effects thus comes from the yet uncalculated radiative 
corrections to the twist 3 contribution: a +30% correction would increase the form factor 
by about 10%. 

Our final result for the spectrum dB(B — ► 7ieu)/dq 2 as a function of q 2 using the QCD 
sum rule motivated pion distribution amplitude and taking central values of the input 
parameters is shown in Fig. 3 (solid curve). We evaluate the light-cone sum rules for 
q 2 < 17GeV 2 , where both the twist-expansion and the contribution of higher states are 
well under control. At higher values of q 2 the decay rate is strongly suppressed by the phase 
space factor; for comparison, we show (dotted line) the spectrum calculated in the vector 
dominance (VDM) approximation in this region, with the coupling gBB*-* = (29 ± 3) ]5|. 
The uncentainties are illustrated by dashed lines. The spectrum decreases monotonically 
with q 2 , which is a consequence of the broad pion distribution amplitude of The 
asymptotic pion distribution produces a different shape, see the dash-dotted curve, with 
a maximum around q 2 ~ 16 — 18 GeV 2 . We repeat that these two distributions present 
two extreme possibilities and the ambiguity will eventually be removed. And vice versa, 
measuring the spectrum in B — ► 7rez/ decays can distinguish between different shapes of 
the pion distribution amplitude. 

To summarize, in this letter we have calculated the radiative correction to the light-cone 
sum rule for the semileptonic B — > ixev form factor. We have studied its behaviour in the 
heavy quark limit, which, as we believe, teaches us that factorization of "hard" and "soft" 
subprocesses is generally valid in the heavy quark limit. From the numerical analysis of 
the sum rule, we obtain the form factor at q 2 = 0: 

/+(0) = 0.25 ±0.03 ±0.01 ±0.01 (20) 

assuming the asymptotic pion distribution amplitude and 

/ + (0) = 0.30 ±0.03 ±0.01 ±0.01 (21) 



using the QCD sum rule motivated distribution [19 1 . The first error comes from the 
variation of the Borel parameter within 5 GeV 2 < M 2 < 8 GeV 2 , the second one from the 
error in the quark condensate and the third one from the combined uncertainties in 
and / 37r . Combining all errors in quadrature and averaging over the two different leading 
twist distributions, we get /+(0) = 0.28 ± 0.05 as our final result. 

We also give a simple parametrization of the q 2 dependence. For the asymptotic distri- 
bution we find: 

2 0.25 ±0.03 

U ^ ' ~ 1 - 1.72 q 2 /m 2 B + 0.716 g 4 /m|j' 



and for the QCD sum rule distribution ||19||: 



2 0.30 ±0.03 

' ~ 1 - 1.32 q 2 /m 2 B + 0.208 q 4 /m%' ( ' 

Both representations reproduce the exact light-cone sum rule results to within 1% for 
< q 2 < 17 GeV 2 . The influence of the other input parameters on the g 2 -dependence is 



negligible, so we do not give errors in the denominators in (|22|) and 
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The total combined uncertainty of ±0.05 can be reduced to ±0.03 by calculating ra- 
diative corrections to the asymptotic two-particle distributions of twist 3 and from more 
detailed information on the pion distribution amplitude when it becomes available. 

Acknowledgements: E.B. received support from CICYT research project AEN95-0815. 
V.B. is grateful to G. Korchemsky for the discussion of infrared factorization in heavy 
quark decays and to the UAB and the CIRIT for financial support during his stay in 
Barcelona. He also thanks the IFAE for kind hospitality. P.B. acknowledges discussions 
with A. Khodjamirian, S. Weinzierl and O. Yakovlev. We thank M. Lavelle for reading the 
manuscript. 



References 

[1] R. Akhoury, G. Sterman and Y.P. Yao, Phys. Rev. D 50 (1994) 358; 
R. Akhoury and I.Z. Rothstein, Phys. Lett. B 337 (1994) 176. 

[2] A. Szczepaniak, E.M. Henley and S.J. Brodsky, Phys. Lett. B 243 (1990) 287; 
G. Burdman and J.P. Donoghue, Phys. Lett. B 270 (1991) 55. 

[3] V.M. Braun, talk given at Brussels EPS HEP (1995); published in Brussels EPS HEP 
(1995) 436 ( Eiep-ph/9510404D . 

[4] V.M. Belyaev, A. Khodjamirian and R. Ruckl, Z. Phys. C 60 (1993) 349; 

A. Khodjamirian and R. Ruckl, Talk presented at 28th International Conference on 
High-energy Physics (ICHEP 96), Warsaw, Poland, July 1996; published in ICHEP 
96, 902 ( |hep-ph/961036?D . 

[5] V.M. Belyaev et al, Phys. Rev. D 51 (1995) 6177. 

[6] A. Khodjamirian et al., Preprint WUE-ITP-97-015 (|hep-ph/ 9706301 ). 

[7] E. Braaten, Phys. Rev. D 28 (1983) 524. 

[8] V.L. Chernyak and A.R. Zhitnitsky, JETP Lett. 25 (1977) 510; Yad. Fiz. 31 (1980) 
1053; 

A.V. Efremov and A.V. Radyushkin, Phys. Lett. B 94 (1980) 245; Teor. Mat. Fiz. 44 
(1980) 157; 

G.P Lepage and S.J. Brodsky, Phys. Lett. B 87 (1979) 359; Phys. Rev. D 22 (1980) 
2157. 

[9] V.L. Chernyak and I.R. Zhitnitsky, Nucl. Phys. B345 (1990) 137. 
[10] A. Ali, V.M. Braun and H. Simma, Z. Phys. C 63 (1994) 437. 
[11] G.P. Korchemsky and G. Sterman, Phys. Lett. B 340 (1994) 96. 
[12] E. Bagan et al., Phys. Lett. B 278 (1992) 457. 
[13] T.M. Aliev and V.L. Eletskii, Yad. Fiz. 38 (1983) 1537. 



10 



[14] J. Flynn, Talk given at Lattice 96, 14th International Symposium on Lattice Field 
Theory, St. Louis, USA, June 1996; published in Nucl. Phys. Proc. Suppl. 53 (1997) 
168 ( |hep-lat/96100T0|) . 

[15] I.V. Musatov and A.V. Radyushkin, Phys. Rev. D 56 (1997) 2713. 

[16] J. Gronberg et al. (CLEO coll.), Preprint CLNS-97-1477 (|hep-ex/9707031|) . 

[17] V.L. Chernyak and I.R. Zhitnitsky, Phys. Rep. 112 (1984) 173. 

[18] S.V. Mikhailov and A.V. Radyushkin, Phys. Rev. D 45 (1992) 1754. 

[19] V.M. Braun and I.E. Filyanov, Z. Phys. C 44 (1989) 157. 

[20] V.M. Braun and I.E. Filyanov, Z. Phys. C 48 (1990) 239. 



11 



